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INTRODUCTION 


The purpose of this document is to provide the Flight 
Test Engineer with a procedural document for the planning 
and conduct of a flight test evaluation of any inertial 
navigation system (INS). Included in this document is a 
description of the types of inertial navigators that will 
be tested at the Flight Test Center in the foreseeable 
future and guidelines on how to plan the testing from an 
initial estimate of the required flying hours to the writing 
eof the final report. The detailed test procedures required 
to collect the necessary data are covered including sample 
flight cards. Data collection and analysis are discussed 
using actual test data to substantiate the techniques 
addressed. Methods of presenting the data in various re- 
ports are discussed. 







The appendixes contain information on the two computer 
programs used to analyze the test data, NAVAN and CEPLOT. 
The program description and users guide is contained in 
appendixes A and B, with check cases in appendix C. 







THEORY OF OPERATION 


‘Inertial space is defined as that space where Newton's 
laws of motion apply. Inertial navigation is based upon 
measurements made with respect to inertial space. An 
inertial systemldetermines the displacement of the carrying 
vehicle from its starting point by measuring the accelerations 
of the vehicle relative to the earth and integrating the 
accelerations with respect to time. 


The basic measuring instrument of an inertial navigation 
system is the accelerometer, an instrument which measures 
acceleration along a single axis. Inside the accelerometer 
is a pendulous mass which is free to rotate about a pivot 
axis in the instrument. There is an electric pickoff which 
converts the rotation of the mass about the pivot axis to an 
output signal. The output signal is fed to a high gain ampli- 
fier and the output of the amplifier is connected back to a 
torquing coil on the accelerometer. When an acceleration is 
present, a current is sent back to the torquer which is pre- 
cisely the amount required to restore the mass to its initial 
position. The accelerometer torquer can be restored with a 
direct current or it can be restored with pulses. Most 
modern systems use pulse restoration because of the ease in 
which the accelerometer output can be processed by a digital 
computer. A simplified block diagram is shown in Figure l. 
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Figure] Acceleremeter 


The current fed to the torquer is proportional to the 
measured acceleration and provides the electrical signal 
which is fed to the navigation computer. The computer 
integrates the acceleration to produce velocity and then 
integrates the velocity to compute distance. If two 
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accelerometers are mounted at right angles to each other 
on a platform which is maintained level with respect to 
the earth and if one of the accelerometers is directed to 
true North, it is possible to determine the distance that 
the platform traveled in the North-South and East-West 
directions. If the platform is "told" where it is initial- 
ly, it can determine where it is on the face of the earth 
at all times. A simple block diagram is shown in Figure 2. 





Figure 2 Accelerometers on Platform 


The accelerometers must be maintained level and in 
proper azimuth~orientation regardless of aircraft attitude 
because even a slight misalignment can introduce serious 
errors. This proper attitude and azimuth orientation is 
maintained by mounting the accelerometers on a gimbaled 
platform with gyroscopes used as sensing elements to con- 
trol platform orientation. A platform which is controlled 
by gyros in this manner is referred to as a stable element. 


A stable element requires three gyros mounted at right 
angles to each other. One gyro senses movement about pitch, 
One about roll, and one about azimuth. Some platforms only 
require two gyros but they are two degree of: freedom types; 
i.e., they are sensitive about two axes. The stable element 
is mounted on gimbals to isolate it from angular motions of 
the aircraft. The operation of the gimbal driving system is 
illustratediim Figure 3. 
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Figure 3 Simptified Single-Axis, Gyre-Stabilized Platform 


The high speed spinning wheel or rotor in Figure 3 resists 
any motion and exhibits an important gyro characteristic call- 
ed rigidity. Any movement about the bearing axis in Figure 3 
generates a pickoff signal which is amplified and drives the 
stable element in an equal and opposite direction to the move- 
ment. 


A practical inertial system requires that the platform be 
stabilized in all three axis of operation in order to retain 
its level orientation regardless of the maneuvers made by the 
aircraft. Figure 4 illustrates a four-gimbal platform configu- 
ration as actually used in an inertial system. The extra roll 
gimbal is provided to prevent the occurrence of a condition 
known as gimbal lock during certain aircraft maneuvers. The 

rn gimbals are oriented so that aircraft attitude and heading 
May be sensed by measuring angles between the gimbals and the 
Platform frame. Synchros transmit this information to the 
attitude indicator and other systems in the aircraft. 
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Figure 4 Feur-Giubal Platform 








The inertial navigation system described up to this 
point is only capable of navigating on a flat non-moving 
earth. This is because the gyros are stablized with 
respect to inertial space and therefore will cause the 
stable element to rotate with respect to the earth as 
the earth rotates on its axis. This is undesirable for 
navigation because the accelerometers will not remain level 
with respect to the direction of gravity. This characteris- 
tic is illustrated in Figure 5. 
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Figure § Space Stabilized Platform 
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When that stable element is mounted in an aircraft the 
same phenomenon occurs at the combined rates of the earth 
rotation plus the aircraft velocity. If the aircraft flies 
North so as to remove the rotation of the earth from the 
sensitive axis of the stable element, the aircraft "sees" 

a continuing pitch maneuver. At the pole, instead of the 
platform being level with the surface of the earth, it would 
now be tilted 90 degrees off level. This characteristic is 
illustrated in Figure 6. 
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Figure § Space Stabilized Platters Mounted in Alreraft 


Obviously this is not a practical system of navigation 
over a spherical, rotating earth so some method of keeping 
the stable element perpendicular to the earth's gravity 
vector is required. By applying a torque to the appropriate 
axis it is possible to maintain the stable element with 
respect to the earth and aligned in agimuth to a known 
reference. Gyro torque is produced by sending a current 
through torquer coils attached to the gyro gimbal. This 
torque causes the gyro to precess at a right angle to the 
applied torque. The computer supplies the current to 











properly torque the gyros. The operation of the platform 
with proper earth rate and vehicle rate gyro torquing is 
= illustrated in Figure 7. 
NORTH. 
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Figwe 7 Earth Stabilized Platform 


The torquing necessary to maintain the stable element 
level with respect to the gravity vector is referred to as 
the transport rate torque. The mechanization of this trans- 
port rate requires a computing loop which is said to be 
Schuler tuned. This computing loop (which is the same for 
. each platform axis, except azimuth) is shown in Figure 8. 
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. Figure 8 Schuler Computer Loop 





The torquing rate necessary to maintain the platform 
level with respect to the surface of the earth is equal to 
aircraft velocity divided by the radius of the earth. The 
system is mechanized so that velocity is multiplied by the 
reciprocal of the earth's radius as shown in Figure 9, 
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Figure 9 Schuler - Tuning 


When an inertial platform is mechanized in this way, 
its response is similar to a pendulum with a length of 3,400 
miles and a period of 84.4 minutes as shown in Figure 10. 
In effect, the center of gravity of the pendulum bob re- 
mains at the center of the earth and the point of suspension 
is at the aircraft. The point of suspension of a pendulum 
can be moved without causing the pendulum to oscillate so 
Schuler tuning allows the platform to be moved about the 
surface of the earth without disturbing the platform level. 
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Figure 18 Schuler Pendulum 
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An important advantage of Schuler tuning is that many of 
the instrument errors are constrained from increasing with 
time and instead are oscillatory in their buildup. An 
example of a constrained acceleration error is illustrated 
in Figure 11. This error is due to an initial level axis 
misalignment of approximately one minute of arc. Note that 
the platform is initially tilted to the left so the accelero- 
meter senses a component of gravity. The system "thinks" it 
is moving and a velocity is developed which torques the plat- 
form in the opposite direction which tends to cancel the 
original tilt error. The system then overshoots level and 
develops an acceleration error in the opposite direction. 
This oscillation will continue until the system is switched 
off or damped by a velocity input from an outside source 
and is referred to as a Schuler oscillation. 
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Figure 11 Constrained Acceleration Errer 


Correcting the platform for transport rate. error allows 
the system to be used for navigation but not very accurately. 
There are three more sources of error that must be corrected 
before the system will navigate accurately. They are 
the Coriolis effect, the oblateness (flattening at the poles) 
of the earth, and centripetal acceleration. 
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The Coriolis effect exists because the earth rotates. If 
the earth were not rotating, a vehicle flying a straight 
ground course from the equator to the North pole would make 
a straight track as seen in Figure 12a. However, because 
the earth is rotating, an observer in space looking down on 
the earth would see that. the airplane really has to fly a 
curved track in space in order to make the desired straight 
ground track over the earth (refer to Figure 12b). Relative 
to space, the airplane must continuously change the magni- 
tude and direction of its tangential velocity. 
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Figure 12 Ceslelis Effect 


Regardless of the direction of the vehicle's horizontal 
velocity, the Coriolis effect appears to be an. acceleration 
to the right in the Northern Hemisphere and to the left in 
the Southern. Although the magnitude of the Coriolis accele- 
ration is small, it can cause a significant navigation error 
when flight times are long. To compensate for the Coriolis 
effect, a correction factor is introduced into the output . 
signals of the accelerometers which is equal to 22Qvsiné, i: 


where 2 is equal to earth rate (15.04 deg/hr), v is velocity er 
(N-S or E-W), and @ is latitude.+ 4 














lperinea in glossary 








ee 





The oblateness of the earth produces a spurious accelera- 
tion in the N-S accelerometer when the plumb line to the 
center of the earth does not exactly coincide with the true. 
vertical as shown in Figure 13. The navigation computer must 
provide an earth radius correction term to the N-S accelero- 
meter to correct for the oblateness. 

\ 
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Figure 18 Spurious Acceleration From Oblateness 


Centripetal acceleration (a,) is a true acceleration with 
respect to inertial space that is created by a vehicle travel- 
ing about a spheroid on any course except a great circle 
course. Figure 14 shows an aircraft flying East along a line 
of constant latitude of radius r'. The center of curvature 
of the path of the aircraft does not pass through the center 
of the earth so a centripetal acceleration is generated. Since 
the path shown is only to the East and is in the Northern 
hemisphere, a South acceleration component is generated. 

To correct for the centripetal acceleration, a correction 
factor equal to -v2/r'tan? must be added to the output 
signal of the North-South accelerometer. 





| 
: 
i 


ee hs alee ial le Aa ae 





~ = 
"ean aae =” 


a, = CENTRIPETAL ACCELERATION a, = CENTRIPETAL ‘ACCELERATION 
SENSED BY ACCELEROMETER 


Figure 14 Contripetal Acceleration 
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So far in this discussion, only the local vertical, North 
pointing inertial navigation system has been considered. 

- This type of INS is referred to as a semi-analytical inertial 
system and is the most common system in use today. The 
platform gimbal structure is simple and the computer mech- 

. anization is easy. The semi~analytical system maintains 
the stable element perpendicular to the earth's gravity 
vector at all times. Accelerometer outputs are con- 
verted to velocity and distance. Velocities are used to 
torque the stable element to maintain the platform normal 
to the earth reference. The stable element is aligned 
in azimuth to an azimuth reference. A simple block diagram 
of a semi-analytical inertial system is shown in Figure 15. 
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Figure 15 Seml-Analytical Inertial System 





oa eins 
Se eS ERE LE ARETE 


iss x issfzdedickcab ainda da =n edt alee 





We hee OS pv 


| 
* 





one? ‘of tit Daath! disadvantages and sources of navigation 
error’ of the semi~analytical inertial platform is the re- 
quirement to torque ‘the gyros to maintain the platform per- 
pendicular to’ the earth's gravity vector. This problem can 
be overcome by ‘using an analytical navigation system which 
does not’ require gyro torquing. Because the analytical 
navigation’ ‘platform remains fixed in space and rotates with 
respect to earth, the accelerometers must sense the gravity 
component and vehicre’ accelerations. For navigation purposes 
only the vehicle accelerations are used and the srayitational 
accelerations must be cancelled out. Calculating the earth's 
gravitational acceleration is an extremely difficult problem 
and requires that an enormous amount of data be stored in 
the computer memory. Acceleration corrections and present 
position information are developed in the computer and coordi- 
mate converter. A simple block diagram of the analytical 
inertial, res is shown in Figure 16. 
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Figwe 1g Analytical inertial System 
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A third type of inertial navigation system is the strap- 
down system. As the name implies, the strap-down system does 
not have a conventional gimbal mounted stable element. In- 
stead, the gyros and accelerometers are mounted directly on 
the vehicle frame. The computer in the strap-down system 
must compute the B matrix computations necessary to specify 
vehicular attitude with respect to an inertial reference 
frame. This computation is usually in direction cosine 
notation, direction cosines being any space vector represent- 4 
ed by three cosines. 


The coordinate converter utilizes inputs from the 
accelerometers and the B matrix to resolve accelerations 
in an inertial reference. A position computer converts 
inertial accelerations and altitude information to 
cartesian coordinates representing vehicle position in 
inertial space. A vector computation then provides 
outputs in latitude and longitude. 


Severe torquing requirements are placed on the 
gyros in a strap-down system because the gyros change 
attitude at the same rate as the vehicle. A simple block 
diagram of a strap-down inertial system is shown in Figure 
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Figure 17 Strap-Down Inertial System 






ERROR ANALYSIS 


GENERAL 





The magnitude of the errors in an inertial system may 
be considered as random variables at any specific flight 
time. In the analysis that follows, it is assumed that all 
errors are independent, linear, and with negligible latitude 
effects on Coriolis and earth rate terms. 

ERROR SOURCES 

Gyro: 
Bias (drift) 
Proportional bias (g sensitive)! 
Anisoelasticity (g? sensitive)! 
Torquing error 
Random drift 
Scale factor 

Accelerometer: 
Bias 
Scale factor 
Cross coupling! 
Vibropendulous! 

Platform: 
Initial level 
Initial azimuth alignment 
Servo errors 
Component non-orthogonality 

Gyro 


Accelerometer 
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Computer 
Round off error (register) 
Truncation error (integration approximation) 
Commutation error (direction cosine) 
Pick off error (acceleration roundoff error) 
General: 
Geophysical data 
Guidance equations 
Target locations 
The local vertical system will be studied in detail be- 
cause it is a conventional system that is relatively linear 
and, consequently, can be handled with Laplace transforma- 
tion. The understanding of the local vertical system is basic 
to all systems operating near the earth. As shown in Figure 
18 and 19, errors can originate from various sources. For 
this discussion only error sources of significant magnitude 


(greater than 2,000 feet/hour) will be considered and all error 
sources will be considered to be step functions. 
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Figure 19 Error Sources in Azimuth Axis 


VELOCITY ERRORS 


A velocity error (€,) will cause the platform to torque 
out of level. The out of level platform will sense a com- 
ponent of gravity which it interprets as an acceleration. This 
signal is integrated into a velocity of opposite polarity to 
the initial velocity error. These velocities oscillate at the 
characteristic Schuler period of 84 minutes. Refer to Figure 
18 where €, is a velocity error inside the Schuler loop. From 
Figure 18, the block diagram shown in Figure 20 can be con- 
structed. 





Figure 20 Velocity Error Loop 
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The transfer function of this error in the system is: 





&qa.™ Ey S2 x2 


converting to time domain 


sindt 


Equation 1 
= velocity error in ft/sec 
Eg = distance error in ft 
g = gravity constant = 32.2 ft/sec? 
r = earth radius = 2.09x107 ft 
t = time in seconds 
Equation 1 can be converted to a more useful form. 
€qg = -l3eysinut ; Equation 2 
€qg = distance error in miles 
w = 4.46 radians/hr 
t = time in hours 
The velocity errors inside the Schuler loop generate 
position errors like the ones shown in Figure 21. velocity 


errors are easily diagnosed if the inertial velocities can 
be compared with actual ground reference velocities. 






€q= 13 €, ‘sin wt 


100 120 


MINUTES 


Figure 21 Effect of Velocity Errer on Position Errer 
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ACCELERATION ERRORS 


Acceleration errors (cq) are integrated into an erroneous 
velocity which, thru the Schuler loop, torques the platform 
out of level, The out cf level accelerometers sense a com- 
ponent of gravity which is opposite in polarity to the accel- 
eration error. For ease of analysis, the input error is assumed 
to be a step function. The vector sum of the acceleration 
error and the gravity error oscillates at the Schuler frequency. 
The error in computed distance is the double integral of this 
acceleration error. In block diagram form, this dynamic csn- 
dition is shown in Figure 22. 





Figure 22 Acceleration Errer Loop 


The transfer function of this error is: 
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Multiplying by > for a step input gives 
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Converting to the time domain 
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Equation 3 can be converted to a more useful forn. 


Eg = 3.45€q(1-cosut) Equation 4 
€q = distance error in miles 
€q = acceleration error in ft/sec? 
w = 4.46 radians/hour 
t = time in hours 


Acceleration errors generate position errors as shown 
in Figure 23. Note that position error does not increase 
with time but is constrained by the mechanization of the 
Schuler loop. Acceleration errors are difficult to isoiate 
in any kind of operating environment other than a laboratory 
because the accelerometer outputs are rarely available 
external to the inertial platform for measurement by the 
instrumentation system. The normal method of. examining the 
accelerometer outputs is to differentiate the velocity 
outputs. 
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Figare 23 Effect ef Acceleration Errer on Position Errer 











LEVEL GYRO DRIFT ERRORS 


Level gyro drift (e€,) is one of the most common sources 
of inertial error. The error is usually the result of an 
improper gyro bias being applied to one or both of the level 
axis gyros. The block diagram for the error is shown in 
Figure 24, 





Figure 24 Level Gyre Drift Erver Loop 


When an error enters the system as a level gyro drift 
error the platform is unleveled. The accelerometers pick 
up a component of gravity which is integrated into velocity. 
Through the Schuler loop the erroneous velocity computation 
torques the level stabilization gyro in opposition to the 
error signal. However, in an undamped pure inertial system, 
the velocity computation overshoots the error signal and 
oscillates as a negative sine wave at the Schuler frequency 
about a ramp function of time. The transfer of the error is: 


ot a drama 
Eg S(S74+12) 
Multiplying by 1/S for a step function input error gives 
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Converted to the time domain 


cq * egr (t-SiRAE, 





Equation 5 


€q = distance error in feet 

€g = gyro drift in radians/sec 

x = earth radius = 2.09 x 10” feet 
t = time in seconds 


Equation 5 can be converted to a more useful form by con- 
verting units. 


Eq = 60eg (t-0.22sinwt) Equation 6 
Ep = gyro drift in deg/hr 
t = time in hours 
w = 4.46 radians/hr 
€q = distance error in miles 


Level axis gyro drift generates position errors as shown 
in Figure 25. 
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Figure 26 Effect ef Level Gyre Drift on Pesition Errer 
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AZIMUTH DRIFT ERRORS 


The most significant cross-coupling error between axis 
is that of azimuth drift rate (5,). This error (Figure 19) 
is integrated by the azimuth stabilization loop resulting 
in an azimuth misalignment angle.t When misaligned, the east 
gyro picks up a component of the earth's rate of rotation 
which torques the stable element out of level. fhe re- 
sulting gravity error sensed by the accelerometer is 
integrated into velocity and torques the platform out of 
level opposite to the earth rate torque through the Schuler 
loop. The computed velocity signal overshoots and oscil- 
lates about the earth's rate torque. The block diagram 
for this error is shown in Figure 26. 
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Figure 26 Azimuth Drift Errer Loop 


From this, the error transfer function is derived. 
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Multiplying by 1/S for a step function error input gives 
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. €qg = sgncosts(=— 


Converting to the time domain 
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_ i-cosirt Equation 7 
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2 a 

= distance error in feet 
azimuth drift rate in radians/sec 
earth rotation rate = 15.04 deg/hr 
local latitude 
earth radius=2.09 x 107 ft 


= time in seconds 


Equation 7 can be reduced to a more useful form. 


Eq = 7.86costs,t* 


Equation 8 


€qg = distance error in miles 


6a-= azimuth gyro drift in deg/hr 


$@ = local latitude 


t = time in hours 
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This error is predominately in the North-South direction 
because the effect of earth rate upon the north gyro will 
be in error by the cosine of the azimuth misalignment. The 
cosine does not change significantly for small error angles. 
The distance resulting from agimuth gyro drift is shown in 
Figure 27. Note that the error is in the North-South axis. 
only and that it increases in magnitude as the square of time. 
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Figure 27 Effect of Azimuth Drift on Pesition Errer 


VELOCITY ERRORS OUTSIDE THE SCHULER LOOP 


Velocity errors entering from outside the Schuler loop 
(€yo in Figure 18) affect the rate of the second integration 
and cause the distance computation error tg to have a straight 
line increase when the error is constant. 


INITIAL AZIMUTH MISALIGNMENT 


Initial azimuth misalignment (e€, in Figure 19) contributes 
to error in resolving horizonal accelerations. It also intro- 
duces errors by causing the east gyro to sense earth rate. 


INITIAL LEVEL MISALIGNMENT 


In addition to accelerometer null uncertainties, there 
may be acceleration errors because of initial level misalign- 
ment errors (€9, €4 in Figure 18) due to a dead band in the 
servo null of the platform gimbals. This error may be written 
as €, = gsinég for a level misalignment in the pitch axis 
where g is the acceleration due to gravity in feet-per-second- 
per-second and €, is the level misalignment of the pitch axis 
in radians. An error in the roll axis would be written as 
Ea = gsineg. 
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TEST PLANNING 
ADVANCE PLANNING 


When the Flight Test Center is assigned as the Responsible 
Test Organization for any flight testing, a certain amount of 
advance planning is required. If the testing is to include 
an inertial system, an estimate of the number of flights and 
flying hours required to evaluate the system will be made. 


The evaluation of the INS should be divided into two 
separate areas, accuracy and operational suitability. To 
predict the accuracy of an INS with an 85 percent confidence 
in the validity of that prediction requires a minimum of eight 
dedicated and valid flights on an instrumented aircraft for 
each alignment mode that has an accuracy specification written 
against it. The eight flights per alignment mode produce a 
Circular Error Probable~(CEP) prediction. Each flight length 
should be in excess of one Schuler period of one hour and 24 
minutes. It is desirable to exceed two Schuler periods but 
that is not a requirement as long as the test flight lengths 
are compatible with the operational mission requirements of 
the aircraft. Operational suitability does not require dedi- 
cated test flights but adds 0.1 flying hours to every test 
flight. 


To determine the inflight accuracy of an inertial system 
with a reasonable confidence level requires an on-board instru- 
mentation system. The inertial parameters must be measured 
with a degree of precision that will allow testing to the 
defined requirements. The specific parameters to be measured 
for each inertial system will change somewhat from aircraft 
to aircraft but the general requirements do not change. The 
following list is the general measurements required for quanti- 
tatively testing an inertial system: 


Inertial longitude 
Inertial latitude 
- Inertial elevation 
Inertial ground speed 
f Inertial ground track 
er Inertial velocity in N-S direction 
Inertial velocity in E-W direction 
Inertial vertical velocity 
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Inertial wander angle 

Inertial roll angle 

Inertial pitch angle 

Inertially computed magnetic variation 
Inertially computed wind speed 

Inertially computed wind direction 

True airspeed input to inertial system 
Magnetic~heading input to inertial system 
Pressure altitude input to inertial system 
Roll angle from auxiliary reference system 
Pitch angle from auxiliary reference system 
Magnetic heading from auxiliary reference system 
Radar range 

Radar bearing 

TACAN range 

TACAN bearing 

Present position update command 

Enter visual fix command 

Aircraft angle of attack 

Aircraft vertical acceleration at the center of gravity 
Tone and event command 


As soon as the location of the test facilities are 
defined at Edwards AFB, the systems test engineer should 
determine if the alignment coordinates and local magnetic 
variation is available for all of the proposed parking spots 
for the test aircraft. The local magnetic variation should 
be remeasured for each of the parking spots approximately 
every four years. This is because the magnetic variation 
in the Edwards AFB area changes at the rate of 1.5 minutes 
per year or 0.1 degrees every four years. 


DETAILED PLANNING 
Test Information Sheet (TIS): 


After the preliminary test planning has been calcu- 
lated, documented, and provided to the project engineer, the 
systems analyst (test engineer responsible for the analysis of 
the flight test data on the inertial system) should begin the 
detailed test planning. This planning is documented on a 
for inclusion in the test plan. The TIS form shown in Figure 
28 has the major topics identified that are to be included in 
the documentation. These topics will be discussed individually 
in the following paragraphs: 
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References 

Test Item Description 
Test Objectives 
Success Criteria 

Data Requirements 
Test Procedures 


Support Requirements 


AFFTC tun 79 261 REPLACES AFFTC FORM 0-128, JUN 73, WHICH WILL BE USED. 


Figure 28 Test Information Sheet 
a 








1.0 References: The first reference to be listed 
is the Air Force management document that is used to generate 
the aircraft or system contract. This document will detail 
the aperational requirements that were intended to be satisfied. 
The next reference document will be the contractors system 
specification document. In this document, the contractor 
tries to quantify the operational requirements contracted 
for by the Air Force. The third reference document 
should be the detailed INS specification document published 
by the manufacturer of the INS. The fourth reference docu- 
ment should be the aircraft flight manual showing how to 
Operate the INS. All of the maintenance documents used at 
the AFFTC to maintain the system should also be referenced. 


2.0 Test Item Description: The description should 
not be of an INS in general but rather, what makes up this 
particular inertial system. Identify all of the aircraft 
components that are considered to be a part of the inertial 
system and will be included in the evaluation. Describe the 
alignment modes and operating modes. Briefly discuss any 
software that will be considered a part of the INS. 


3.0 Test Objective: The test objective will be 
much the same for any INS evaluated. It will normally 
be "To determine the operational accuracy and usability of 
the Inertial Navigation System." 


4.0 Success Criteria: The accuracy portion of the 
testing will be complete when sufficient test points are 
available to give a reasonable estimate of the INS accuracy. 
This estimate requires a minimum of eight valid-data flights 
for each mode of operation and/or alignment. If any 
development changes are made to the hardware or software 
after the accuracy testing has been accomplished and if those 
changes could affect system accuracy, then sufficient re- 
testing must be accomplished to demonstrate the impact of 
the change. Operational suitability testing will continue 
on a ride-along basis until the end of the test program. 


5.0 Data Requirements: The end point data for all 
valid data flights is normally presented on a circular error 
plot with the vertical axis showing North-South error and 
the horizontal axis showing East-West error. The accuracy 
requirement is shown as a circle whose radius is equal to 
the specification CEP. Each flight terminal error point is 
normalized to the time duration of the specification and 
plotted. The average value of all the plotted points is 
calculated and plotted. The predicted CEP based on the test 
data is derived and plotted according to a formula that is 
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used by the AFFTC to evaluate all inertial systems tested 

at the Center. If the contractor uses a different technique 
for calculating CEP, that value should also be calculated and 
plotted. 


If in-flight accuracy data are obtained, the data will 
be presented as an error time-history plot for each flight. 
Multiple flights in an alignment mode will be used to generate 
a CEP time-history plot. A separate CEP plot will be made for 
each alignment mode evaluated. An 85 percent confidence limit 
should be shown on each CEP time-history plot along with the 
specification for the INS. 


Ground test data should be obtained and compared with the 
specifications called out in the maintenance technical publi- 
cations for ground checking the platform. The relationship 
of ground accuracy and in-flight accuracy should be docu- 
mented. 


6.0 Test Procedures: This section will contain the 
overall test plans for the inertial platform. It should be 
as detailed as possible at the time the TIS is prepared and 
should start with ground testing and conclude with in-flight 
testing. Specific test techniques will be discussed for 
various types of testing later in this report. 


7.0 Support Requirements: This section should detail 
all of the services furnished by organizations other than 
Systems Engineering. Included should be a list of the instru- 
mentation parameters, range support, data support, radar track- 
ing support, and photographic support. 


Data Collection: 


Collection of the INS data is the responsibility of 
the systems test engineer assigned to evaluate the INS system. 
The test engineer must plan in great detail how the test 
points are to be flown. After he has planned how to obtain ) 
the test points on each flight, he should make up the test 
cards for this part of the testing. He should discuss 
the test cards with the Project Engineer to make sure that 
all test points and cards are practical. A test support 
summary should be written up for each test flight detailing 
flight speed and length, photo and/or safety chase require- 
ments, radar tracking requirements, photo theodolite tracking 
requirements and instrumentation support requirements. This 
planning must be complete before scheduling a test mission 
(approximately two weeks before flight). 











Instrumentation support requirements planning should 
produce a prioritized list of parameters which must be opera- 
ting for data analysis on each flight. The list should in- 
clude measurements of aircraft configurations, aircraft 
attitude, and aircraft dynamic conditions so that it can be . 
determined what the dynamic conditions of flight were at 
specific times during the flight. The planning should also 
identify any special instrumentation pre-flight and post 
flight requirements. 


A complete history will be kept on each inertial 
platform in the test program. When an inertial platform 
is initially assigned to the test program, a platform history 
log like the one shown in Figure 29 is initiated. Every 
hour of platform operation during the test program should 
be documented, regardless of the reason for the operation. 
This requires the cooperation of the ground maintenance 
personnel. Control procedures must be worked out with the 
maintenance personnel and the project engineer on each 
test program to make sure that none of the data is lost. 


The accuracy of the INS platform boresight should be 
determined and documented at some time during the test program. 
This is necessary to prove that INS data obtained on any parti- 
cular aircraft is valid and does represent what can be expected 
out of the using command aircraft providing they are also within 
tolerance on the boresight alignment. The boresight check should 
be made as early as possible in the test cycle. If the plat- 
form is out of tolerance, all data obtained prior to correcting 
the boresight is invalid and can not be used to demonstrate : 
the system capability. 


Inertial data are obtained during all-weather testing 
and should be analyzed. However, that test requirement should 
be considered as being in addition to the baseline data that 
will be obtained at Edwards AFR. The reason for this is that 
if there are any observable differences between the data be- 
cause of the difference in environment, there must be enough 
data at each condition to prove and present the difference. 


Keep in mind that the all-weather test aircraft will be tested 
at a remote site so it may be impossible to analyze the data : 
in a very timely manner. In fact, there is a good probability 
that the test engineers that accompany the all-weather test 

aircraft to a remote test site will not have the background 

to analyze the INS test data. Their primary responsibility id 
is to conduct the test, not to analyze subsystem performance. 
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Figure 29 Inertial Platform History Log 
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Ground Testing: 


Sufficient ground testing should be accomplished to 
establish the relationship of the ground accuracy to the 
in-flight accuracy and to validate the ground maintenance 
procedures called out in the T.O.'s. This testing should be 
accomplished by the engineer responsible for analyzing the 
INS performance but could be accomplished by others. As an 
example, maintenance will often perform a ground drift run to 
check out the system. The engineer should obtain and use the 
data from the drift run. 


The length of a ground drift run should never be 
less than 84 minutes. A minimum of eight ground drift runs 
should be made for each alignment mode during the course of 
the test program to obtain a confidence level of at least 
85 percent in the validity of the data. If in-flight 
accuracy data are obtained, it is a good policy to obtain 
drift run data before and after the flight using the same 
alignment mode as for the flight. However, the platform 
should be allowed to cool down to ambient temperature 
between shutdown and the next alignment and that requires 
a minimum of 12 hours. It is also preferable to start the 
alignment on a different heading than the shutdown heading. 
This will insure that the platform will have a reasonable 
amount of shutdown error to correct during alignment. Align- 
ing on various headings helps to identify the effectiveness 
of the heading bias correction constant. A sample INS 
ground drift run card is shown in Figure 30. 


Operational Testing: 


For operational testing, the data that will he plotted 
is normalized endpoint data. The inertial platform 
history log shown in Figure 29 has all of the data required 
to produce endpoint data plots. The flight data on that log 
will be obtained by the systems test engineer from the : 
flight crew at postflight debrief. Sample alignment and shut- 
down cards are shown in Figures 31 and 32. These cards should 
be prepared by the systems test engineer and given to the 2 
project engineer at least one day prior to flight. The a 
REMARKS entry shown on Figure 31 is to document the wind con- 
ditions during the alignment, any observable vibration or 
movement of the aircraft during the alignment, and whether the 
power source was switched from ground power to aircraft power 
during the alignment. 
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INS Ground Drift Run Test Cards (Cards 1 and 2) 


Figure 30 


L 


apnzF8uo7 | apnazie7] | aut, pesdety 








Z@ auvo - NNW Ld41ud GNNOUD SNI 








our 


Jjo wa3sks [eTIAeuT 3sy. uaANy 
Xy 


dds aNd 
AQPIOTAA wWeAsks [eTIAaUT pz0d9y 
epnqTguoy apnqti.ey 
Sa}eUTPI0OD weIzSks TeTIIsUT p1rodsy 
~~~ Bu yreeg ———Suy 


sa JeutTpioos 
uaop3nys 03 3uz1eeq pue e3ue1 paznduos wezsks 


YUVA OWN 
UOTIETFACA JypJeuseu weysks 
3PH 


Suypecoy we 3sts 


@pnz ;38u07 epnz tie Ty 


SaeUTpPIOOD UAMOP NYS 


NMOGLAHS IVILYANI 


Inertial Shutdown Test Card 


Figure 32 
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Inertial Alignment Test Card 


Figure 31 











The flight data will be presented on data plots that 
reflect the alignment mode. To demonstrate any type of rapid 
alignment, the INS should be allowed to soak at ambient con- 
ditions for a minimum of 12 hours. before starting the align- 
ment. If the ambient soak conditions can't be met, then the 
platform should be aligned in the standard alignment mode 
(gyro-compass for a semi-analytical system). 


If the inertial platform has to be shut down in 
flight for any reason, the data recorded before platform 
shutdown can not be used for endpoint accuracy data. How- 
ever, if the platform is brought back up with an in-flight 
alignment, that data can be used as end-point data for the 
in-flight alignment accuracy plot providing it is possible 
to define the accuracy of the aircraft position versus the 
entered initialization ccordinates. 


Often during a test program it will be necessary to 
enter a present position correction to the INS during flight. 
For a non-instrumented aircraft, that almost always prevents 
the test engineer from using the data from that flight for any- 
thing other than a qualitative description of how well it 
worked. This is because it is necessary to know where the 
aircraft was immediately before update, how much error was in 
the INS at that time, and how much of the error was removed 
by the update in order to use the data from the flight in 
any kind of meaningful analysis. All of this information 
must be time correlated to the alignment and shutdown data. 


Many inertial systems tested in recent years are 
bounded-error systems. This means that position and velocity 
parameters from the INS are mixed with data from other sources 
by a statistical filter to arrive at a computer estimate of 
position and velocity. This computer estimate is then fed 
back to the INS to try to keep the difference between them as 
small as possible. It is not really possible to isolate the 
INS accuracy because the accuracy of the other systems show 
up in the results. 


Air Force accuracy requirements are nearly always 
specified for operational conditions. However, the con- 
tractor normally demonstrates specification compliance under 
a limited set of carefully controlled conditions. The 
initial responsibility for quantitative testing under 
operational conditions falls to the Flight Test Center and 
the systems test engineer should lay out his testing require- 
ments with that in mind. Alignments should be made in all 
available modes and headings. Maneuvers should be specified 
on a significant number of flights that will stress the inertial 
Platform with the same type loads expected operationally. 

The systems test engineer must personally prepare the 





alignment and shutdown cards for each flight and review 
all flight cards to make sure that the flight will represent 
a fair operational challenge to the INS. 


Inflight Accuracy Testing: 


The data objective of inflight accuracy testing is 5 
to produce a single plot that compares the predicted CEP 
based on actual test data with the specification CEP for 
each alignment mode evaluated. It is desired to have a 
Minimum of eight flights of valid data for each mode to 
calculate a CEP with a reasonable confidence level. 


For inflight accuracy testing of current inertial 
systems, an instrumentation system is usually required and 
is always highly desirable. The system should record the 
on-board position coordinates with respect to a time that 
is correlateable to ground tracking data. As inertial 
accuracy flights are dedicated flights, the instrumentation 
recorder should be turned on and off at the desired test 
point intervals. An instrumentation event tone should be 
transmitted briefly during each recorder-on cycle so that 
each segment of recorded data can be individually checked 
for time correlation. 


Range tracking accuracy is an important considera- 
tion for inertial accuracy testing. For systems that have 
been tested at Edwards in the past, the accuracy of the space ° 
positioning radar has been sufficient if the radar grazing 
angle is kept a minimum of two degrees above the horizontal 
plane of the radar antenna. The azimuth angle accuracy of 
the radar is specified as 0.2 mils which is 1.2 feet per mile 
of range. Figure 33 is a plot of the radar azimuth accuracy 
as a function of range from the tracking antenna. The minimum 
MSL altitude of the aircraft required to obtain this accuracy 
is shown at 25 mile increments. The ranging accuracy for 
the radar is specified as $12 feet at all ranges. 


The flight profile for an accuracy flight can be very 
flexible but because of the amount of data required, it is 
recommended that a rectangular pattern be flown along the 
lines of constant latitude and longitude with the length 
of the legs being as long as possible without exceeding the 
altitude and accuracy restrictions. High rate turns and pitch 
maneuvers should be programmed into the individual legs to 
task the inertial system as much as possible. The systems 
test engineer should monitor the test flight from whatever 
location that is available to him that provides him the 


greatest insight to the data being produced. > 














Figure 33 Space Positioning Reader Asiauth Accuracy 
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The data to be collected during the alignment before 
an inertial accuracy flight is shown in Figure 34. Observe 
that the data card requires the crewmember to duplicate many 
of the parameters being recorded by the instrumentation system. 
This is so that if the instrumentation system should fail, 
there would be sufficient data available to keep the Inertial 
Platform History Log current and use for end-point accuracy 
data. Initial alignment headings should be varied so that a 
minimum of two alignments will be performed in each compass 
quadrant. 


The next data point is on the taxi-way just before 
the takeoff end of the runway. There is an inertial checkpoint 
marked on the taxi-way at both ends and the center of the 
runway at Edwards AFB. Photographs of the checkpoints at 
the ends of the runway are shown in Figures 35 and 36. When 
the nosewheel is over the INS checkpoint, the pilot should stop 
the aircraft, place an event mark on the instrumentation tape, 
and write down the time and coordinates. A suggested flight 
card format is shown in Figure 37. 


“ 

During the flight, test points should be recorded 
every five minutes. Because of the number of test points in- 
volved, the number of test cards for inflight accuracy data 
could become unwieldy for an extended mission. It is important 
to lay out the test in such a way as to assure being able to 
correlate data. A sample test card is shown in Figure 38. 

As stated earlier in this report, the flight should last a 
minimum of one full Schuler cycle of 84 minutes so there 
should be a minimum of 17 in-flight test points five minutes 
apart. 


The Inertial Platform History Log should be filled out 
by the systems test engineer from data obtained in the post- 
flight debriefing just as if it was a non-instrumented flight. 
The fact that in-flight accuracy data exists for a flight should 
be noted in the Remarks column. 


REVIEW OF EXISTING DATA 


Normally, the Flight Test Center will not test an INS 
that does not have a considerable amount of test data already 
available on it. That test data may be from the INS manufac- 
turers development testing, from testing conducted by the 
Central Inertial Guidance Test Facility (CIGTF) at Holloman 
AFB, NM, from the avionics development and integration test 
program of the airframe contractor, or from test reports on 
other aircraft with the same inertial system installed. The 
systems test engineer should attempt to obtain as much of this 


INERTIAL ALIGNMENT 


Align mode 

Instrumentation recorder ON 

Function switch to align 

Enter coordinates and reset 

Heat light out 

Record outside air temperature 
Temp 


Alignment complete 


System heading 


HDG 
System magnetic variation 
MAG VAR 
Select NAV position 


- Instrumentation recorder OFF 





Figure 34 Alignment Test Card - Instrumented Aircraft 
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Runway 22 at Edwards AFB 


Figure 35 Inertial Accuracy Checkpoint No. 1 
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Pre-Takeoff Test Card - Instrumented Aircraft 


Figure 37 








data as possible prior to the start of the test program and 
be prepared to compare it with Flight Test Center data as 
that becomes available. 


If the data from all sources does not correlate, the 
reason must be determined. If the data does correlate, it 
might be possible to reduce the number of data flights re- 
quired at the Flight Test Center. However, this is a function 
of the confidence that the systems test engineer has in the 
data, so reducing the number of test flights should not be 
included in any of the planning. 


REVIEW OF AGE AND TECHNICAL PUBLICATIONS 


It is the responsibility of the systems test engineer to 
review all applicable technical publications and evaluate 
the contractor furnished support hardware to be used by 
Air Force maintenance personnel. The technical publica- 
tions should give clear and accurate directions on opera- 
tion, troubleshooting and maintenance procedures for the 
INS. The test equipment should provide assistance in the 
trouble-shooting and repair of the system. 


An adequate evaluation of the AGE and technical publi- 
cations requires the assistance and cooperation of the 
maintenance personnel assigned to the test program. The 
systems test engineer should try to keep the maintenance 
personnel fully aware of what he is trying to accomplish 
and the results from data already gathered. This includes 
showing them the data, showing how it correlates, and 
discretely pointing out the mistakes that maintenance 
personnel may make that can cause the data to be invalid. 
Any recommended changes in the maintenance procedures should 
be discussed with the maintenance personnel before submitting 
to the SPO. 


REPORTS 


The systems test engineer should plan for the reporting 
phase of the INS testing just as carefully as he does the 
rest of the test program. The reports must present all 
knowledge obtained as a result of the testing. 


The first type of reporting that may be required of the 
test engineer is for Air Force Preliminary Evaluations. This 
type of testing is normally a quick overview where the entire 
aircraft is evaluated on a very small number of missions. 

If an INS had gross problems that would make it operationally 
unacceptable, this type of testing should identify it. How- 








ever, the likelihood of detecting subtle design problems or 
accuracy problems is not very great. The data should be 
presented as a written summary of what was observed along: 
with conclusions about the system operation. 


During the conduct of a flight test program, the Deficiency 
Report (DR) is used to document and track deficiencies in the 
inertial system. The specific directions for using the DR 
form can be found in AFFTCR 80-2. Generally speaking, the 
DR is used to identify suspected problems to the SPO and the 
contractor. The systems test engineer will have the 
responsibility of keeping track of the DR's on the INS 
and whatever corrective actions that take place. A DR form 
is shown in Figure 39. 


The progress report is a periodic management reporting 
tool to keep the report addressees informed on the progress of 
the testing and to provide an insight into any problems that 
exist. The progress report should summarize the INS testing 
for the report period, summarize the results, relate that 
to what has been conducted previously, and describe what 
remains to be accomplished. A summary of the deficiency 
reports submitted during the reporting period should be in- 
cluded. If several flights are involved, an INS flight log 
like the one shown in Figure 40 should be included. Note that 
all of the information presented on the flight log can be 
obtained from the Inertial Platform History Log (Figure 29). 


The final report should cover everything that is known 
by the Test Center about the INS. The body of the report 
should present summarized data to substantiate conclusions 
about the operation of the INS. Individual flight data should 
be in a data appendix to the report unless specific data are 
required to substantiate a conclusion made in the body of the 
report. : 
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Figure 39 AFFTC Deficiency Report Form 
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DATA ANALYSIS 


TERMINAL ERROR PLOT 


The terminal error plot shows the performance of an 
inertial system during the test program in terms of shut- 
down error for all flights. A sample plot is shown in Figure 
41. Each shutdown error obtained during the test program is 
normalized to the time specified in the INS specification 
and plotted. The CEP is calculated and plotted as a circle 
of equal probability for each alignment mode evaluated. The 
title block identifies the alignment mode and the time 
period represented on the plot. The specification CEP is 
shown as a circle of radius equal to the specified value. 


The data for the terminal error plot is obtained from 
the flight log (Figure 40). The normalized values of lati- 
tude error, longitude error and radial’ error are calculated 
by the systems test engineer and stored on a computer data 
card. If the test engineer judges that the data from a 
particular flight is invalid, it should be documented in the 
inertial platform history log and flight log and not punched 
onto data cards. 


When a reasonable amount of terminal error data is avail- 
able, the data is run through a CEP calculation. After the 
calculation is complete, the computer examines the distribu- 
tion of the terminal error points with respect to the value 
of the calculated CEP. If the terminal error from any of the 
flights exceed a three-sigma value of the CEP, that flight 
terminal error is temporarly suppressed from the calculation 
and the calculation is done over. The test engineer should 
include all valid terminal error data that is available in 
the CEP calculations. New flight data should be incorporated 
as the inertial evaluation portion of the testing progresses. 


POSITION ERROR PLOT 


If inflight accuracy data are available, position error 
plots showing the INS accuracy as a function of elapsed time 
should be produced. Figure 42 shows a sample position error 
plot produced from actual test data. The plot displays lati- 
tude error, longitude error, and radial error as a function 
of elapsed time. On the plot, north latitude error and east 
longitude error are defined as positive. Radial error is 
always positive. 


The AFFTC computer program that produces the position error 
plot is called the Navigation Analysis Program (NAVAN) ¢ In- 
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Terminal Errer Plot 


Figure 41 
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Position Error Plot 


Figure 42 












puts to the computer for a particular flight are a radar track- 
ing tape of that flight from the Space Positioning Branch, an 
engineering units tape produced from the aircraft instrumenta- 
tion raw data tape for the same flight, and card entered data 
for information that can't be obtained from the other two 
sources. The program merges the data, edits the data accord- 
ing to given editing specifications, writes a data file of the 
selected and edited parameters, performs computations, and 
produces plots of the computations. 


The first step in producing position error data is to make 
a time correlation table for the flight that is to be plotted. 
This is done by comparing the time of the tone cutoff on the 
Visacorder data recorded at the Space Positioning facility 
during flight with the time of the end of the event pulse 
from the onboard data for each test point. If the time 
difference between the two data sources shifts significantly 
(greater than ten percent of the specification drift when con- 
verted to distance) between the first and last tone breaks, 
the data will be unusable unless each tone-break is correlated 
individually. Both data sources are normalized to make time 
equal zero be the time that NAV was selected on the INS control 
panel. 


Next, the data is run through a smoothing routine and 
interpolated to make the data samples from the two sources 
occur at the same time. The smoothed and interpolated data 
is then compared to obtain the latitude error and longitude 
error. Radial error is calculated from the latitude and 
longitude errors. A listing and the plotting tape required 
to produce the position error plot shown in Figure 42 is then 
created. 


CEP PLOT FROM POSITION ERROR 


A single CEP plot that summarizes the results of all of 
the position error plots for each alignment mode is shown in 
Figure 43. The data used to generate the CEP curve are 
obtained from all of the radial error information from all 
of the individual accuracy flights. The 85 percent confi- 
dence limits shown in Figure 45 make the prediction that 
85 percent of all future navigation errors will fall be- 
tween the two curves. The specification line is the accuracy 
that the contractor was required to meet. This number is 
usually a straight-line specification for a finite time-limit. 

Three nautical miles per hour for the first three hours of 
operation is the specification that is illustrated. 
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Figure 43 CEP Plot from Position Error 
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NAVAN Computer Program 
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INTRODUCTION/PROGRAM DESCRIPTION 


The Navigation Analysis Program (NAVAN) is a software package 
developed to provide the test engineer with a method for determining 
the performance of inertial navigation systems. NAVAN was designed to 
analyze INS data obtained from flight test. NAVAN can accept data 
from Radar and ADAS tapes and/or punched cards. After at least three 
flights are merged, the mean, median, R50, R90, CEP and confidence limits 
are calculated and data is plotted. 


Data editing of the Radar data is done by checking radar velocity 
against the maximum value selected by the engineer (default = 900 ft/sec). 
If the velocity exceeds this value the data associated with this time 
is replaced with the data at the last inbound time. Replaced data is 
noted on the listing. 


To match Radar and ADAS data at the exact same times, a linear 
interpolation is made of the ADAS data. 


Data from the different sources will not occur at the same sample 
rate. Because a uniform rate is required for statistical analysis, a 
least squares fit is made of the data and values from all sources are 
established at the same sample rate. 


The initial flight's data are placed on the New Flight File. If 
not the initial flight, the data are added to previous data and are 
placed on the New History Flight File (Old History Flight File on next 
flight). While a time history plot and listing of the data are avail- 
able, it is recommended that at least three flights be obtained before 
a statistical analysis is made. 


THEORY 

The methods to calculate the 50th and 90th percentile and confidence 
limits on the 50th percentile used in the analysis program are presented 
in this section. Two methods are incorporated in the program for 
calculating the 50th percentile value. 

The following definitions are used throughout this section: 

x = latitude error (nm) 

yY ™ longitude error (nm) 


radial error (7x7 + y7) (nm) 


m = number of tests 
i subscript = test, i= 1 thru m 
q = a variable 


fq, = sum of a qy for i = 1 thru m 
ie 









and R90) is based on the Air Standard 53/11.B, 15 August 1968, The 
Specification and Evaluation of the Accuracy of Inertial Navigation 


System. The computations for R50 and R90 are made at each time point in 
the following manner: 


from the m corresponding flights in the sample. 


(CEP only in this program) is based on a paper by L. L. Rosen and D. L. 
Harmer titled "Inertial System Performance Evaluation" which was 
presented at the Third Inertial Guidance Test Symposium at Holloman AFB, 
New Mexico, in 1966. A procedure for calculating confidence limits on 
this CEP is also given in this paper. 


Wg = mean or expected value of q 

074 = variance of q 

G = sample mean of q = Zqj/m 

S*,q = sample variance of q = {Z(q; - q)?}/m 


The first method used to calculate percentiles of radial error (R50 


Suppose, at some point in time, there are wt radial errors (rj, 1, m) 


1. Calculate the geometric me an+ (GM) of the radial errors: 
m 
GM = ViIry where: I = the product 
2. Calculate the root mean square (RMS) of the radial errors: 
RMS = v¥ir*;/m 
3. Calculate the ratio = GM/RMS 


4. Define RATIO = GM/RMS and calculate the R50 and R90 values from: 


R50 = RMS (.7 RATIO + .3) for RATIO < .6 

R90 = RMS (1 + YI —- RATIO) 

R50 = RMS (.7 RATIO + .4 YRATIO) for RATIO > .6 7 
R90 = RMS {RATIO + 1.6 (1 - RATIO?) } iz 


The second method used to calculate the percentile of radial error 


At each time point the percentiles of radial error are calculated | 


= 
R, Sy Va (2507 + Ver” 


where: 


the pth percentile of radial error and Zp = the pth percentile 


pote 8? a zero mean normal distribution. 


Gyr ar Sze Uz are calculated from the following set of formulas: 
Sy = Sy¥m/(m- 1) 7 oy ® Sy¥m/(m - 1) 


Wy =X Wy *¥ 
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n = K?(2 - K). + 1 + (2/ay*) (d* > uy2K* - Hy?) 


A = K?(K? - 1) + (1/oy?) (2uy?K? 2uy? - a?) 
a=ntdA3; b= A/a 
Uz = 1 - (2/9) (1 + b)/a - (40/81) (b2/a2) 
o2 = ¥(279) (1 + b)/a + (16/27) (b*/a*) 
The 50th percentile (CEP) is calculated from: 
Rp = CEP = oyVau*, 
since: 
25 =o for p = 50 


At each time point the 100 (1 - a)* percent confidence limits on CEP 
(L1<CEP<L9) are approximated by substituting the upper and lower 
confidence limits of the means and sigmas calculated below into the 
formulas above: 


The confidence limits on the means are calculated using a “t test": 
|x] - {(tgq/2) (ay/v¥m)} < wy < [x] + ((tq/2) (o4/%m) } 
I¥| - {(tg/2) (oy/v¥m)} < wy < 1K] + ((tq/2) (oy/¥m) } 


The confidence limits on the sigmas are calculated using a "chi 
test". 


{Ss .vn7(m = TV }/ (x, - a/2) < o, < {s, vn7(m = 1) }/x,/2 
{S,¥m7(m = TV }/ (xz - a/2) < oy < {s/m/7(m = 1) }/x,/2 
where: 


ty/2 is the upper 100 (a/2) percent point of the "t distribution" 
with m- 1 degrees of freedom. YX/2 and xX .7 a/2 are the lower and the 
upper 100 (a/2) percent points of the “chi’distribution" with m- 1 
degrees of freedom. 


The lower limit of the CEP (L}) is calculated by using the lower 
limits of mean and sigma to calculate L) = OyvVay,". Should the lower 
limits of the means be negative, a zero is substituted for this lower 
limit in the computations of L). 


The upper limit of the CEP (L2) is calculated by using the upper 
limits of mean and sigma to calculate L2 = dyvay,2°. 


PREPARATION FOR USE 


The program deck of NAVAN will be permanently stored on magnetic 
tape number 04295 at the AFFTC tape library. Prior to a first use, the 
compiled file must be copied to disc and stored as a permanent file as 
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shown in figure Al. Once the program is stored as a permanent file, it 
may be attached and executed. 


SCOPE OF PROGRAM 


This program uses control information from data cards and extracts 
selected test data information from one or certain combinations of the 
following inputs: (1) An Automatic Data Acquisition System (ADAS) tape 
(produced at the AFFTC), (2) A radar tape (produced at the AFFTC), (3) 
Card data containing navigation system time, latitude, and longitude, and 
(4) Card system data containing time, latitude, and longitude position 
errer. The program then merges the data, edits wild points, performs 
various calculations, produces plots on the Calcomp plotter, and produces 
printed output of the computations. 


Generalized flowcharts of the three data sources (ADAS - radar, radar 
card, and system card) are shown in figures A2, A4, and A6é respectively. 
Detailed flowcharts are shown in figure A29. The subroutines, ADADATA, 
RADRD, and CRDDATA, select the times and parameters from the data sources. 
Multiple events on the ADAS tape are read with one start and stop time. 

A wild point check is made on the radar data only, \»ased on a maximum 
aircraft horizontal velocity of 900 feet per second. (Values replaced 
are noted on output listing.) This value may be changed by the user if 
required. Any radar data input at the time that the velocity exceeds 
900 feet per second is replaced with the last inbound value. The radar 
data is then merged with the ADAS or system card data, based on the time 
of the radar data point. A linear interpolation is used to obtain the 
data value at the time of the radar data point. The radar latitude and 
longitude is subtracted from the system latitude and longitude. The 
resulting error is converted to nautical miles and the radial error is 
computed for each data point. A least squares curve fit is then made 
of the data in order to obtain a constant sample rate. The sample rate 
output is a user defined option. 


Note: For valid statistical data, the sample rate should be the same 
for all flights, regardless of the data source (ADAS - radar, radar card, 
or system card). This data is now placed on the New Flight File (tape 16). 
Tape 16 may be optionally printed or plotted. Tape 16 is normally the 
input to NAVMR. This is combined with the file containing data from the 
other flights to be analyzed (tape 17). This merged data is placed on 
the New History Flight File (tape 11), and may also be optionally 
printed. Tape 11 is the input to NAVAN which performs the statistical 
calculations of the data to be analyzed. An important user defined 
variable in NAVAN is the time span over which each calculation occurs. 
This time span should be equal to the time between samples. This is done 
to obtain only one data value per flight number per time span, in order 
to maintain statistical validity for the calculated confidence intervals. 
NAVAN generates both Calcomp plots of the data and a computer listing. 





USER'S GUIDE 


Listed below are the necessary tapes for options being used: 


LIST OF FILES 
DESCRIPTION 


Search file (Internal Operating File) 

Input (Card Input) 

Output (Printer Output) 

Input (ADAS Tape) 

Input (Radar Tape) 

I/O (ADAS Selected Parameters) 

Output (Optional, from NAVAN for Debugging) 
New History Flight File (NHF, New Data Base) 
ae (Statistical Data for Plot, Internal Working 
File) 

Calcomp Plot Tape 

I/O (New Flight File, NFF) 

Input (Old History File, OHF, Old Data Base) 
I/O (Radar Selected Parameters) 

Internal Working File 


The output data from NAVAN consists of an output file computer 
listing, as in figure A8, and Calcomp plots, as in figures A9, Al0, and 
All. The New History Flight File (file 11) contains all of the previously 
processed flights sequentially arranged in order of increasing time. The 
output listing contains a copy of all of the input cards and a copy of 
all of the processed data. This can be used to check the validity of the 
data prior to plotting. 


INPUT CARDS FOR NAVAN 


The cards for the three different input options are described 
separately. 


ADAS - RADAR COMBINATION DATA 


NAVAN gata card #1, figure Al2, is used to pick the input option, 
INOP(2)=#1, for ADAS and Radar tapes, set the slope of the CEP line to be 
drawn on the Calcomp plots and select other program options. The 
NAMELIST format is used to read in data. A card must start in column 2 
with $DATA and be terminated with a $ sign. Three additional input 
options may be specified on this card. Normally 900 ft/sec is used 
internally for wild point checking, but this may be reset to another 
value by placing N900=XxXxx.xx (the new value) on the card. (Care should 
be taken that the actual aircraft velocity does not exceed 900 ft/sec.) 


NAVAN data card #2, figure Al3, contains the number of parameters to 
be selected from the ADAS tape and the parameter ID codes. Note: One 
data word may be made-up of two words containing the most and least 
significant data bits. 


NAVAN data card #3, figure Al4, contains data to correct the time on 
the ADAS tape to the time on the radar tape. This card also contains the 
navigation system start time. 
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NAVAN data card #4, figure Al5, contains the start and stop search 
times in total seconds. 


NAVAN data card #5, figure Al6é, is in NAMELIST format and starts with 
a $NAM1 in column 2 and is terminated with a $ sign. The next entry is 
the start time, in seconds, on the New History Flight File, tape 11, where 
data calculations start. Next, the stop time, in seconds, on tape 11 is 
listed. The time span for calculations is specified and should be 
selected so as to have one data point per time span per flight for the 
best statistical validity (standard time 300 seconds). 


NAVAN data cards #6, #6A, #7, #8, and #9, figures Al7, Al18, A19, and 
A20, are only required for Calcomp plots. Data card #6 has the X-axis and 
Y-axis scale factors for each plot, the initial value for the X-axis, the 
number of plots, and a plot heading code. The number of plots is equal to 
three times the number of runs. The plot heading code must be equal to 1 
for a heading and 0 or blank for no heading. Data card #6A contains the 
heading. Data card #7 has the axes lengths for Calcomp plot #1, figure A9, 
and the starting value for the Y-axis. Data card #8 is the same as data 
card #7, but contains information for Calcomp plot #2, figure Al0. Data 
card #9 is the same as data cards #7 and #8, but contains information for 
Calcomp plot #3, figure All. 


If there is to be more than one data run, cards with the new values 
should be repeated. The end of the data request should be indicated by : 
a blank data card, figure A211. 


SYSTEM CARD DATA 
This is the case where system error data will be entered on cards. “ 


NAVAN data card #1, figure A22, starts with a $DATA in column 2 and 
terminates with a $ sign. The NAMELIST format is used to read in data. 
This card should be coded with INOP(3)=1 and INOP(10)=1, to indicate that 
the data cards will have latitude and longitude in degrees. Other input 
options are the same as for ADAS - Radar Combination Data, figure Al2. 


NAVAN data card #2, figure A23, will be the first data card. One 
card will be read for each data point, and the final data value will be 
~100000.00, figure A24. The following cards after the final data value 
will be the same as for ADAS = Radar Combination Data. 


RADAR CARD DATA 


This is the case where radar tracking data are available and 
system data are on cards. 


NAVAN data card #1, figure A25, starts with a $DATA in column 2 and 
terminates with a $ sign. The NAMELIST format is used to read in data. 
This card will contain the codes INOP(4)=1 and INOP(10)=1, with the other 
input options the same as for ADAS - Radar Combination Data, Figure Al2. 


NAVAN data card #2, figure A26, contains the start and stop times 
for the radar tape, and the navigation system start time. 


NAVAN data card #3, figure A27, is the first data card and should be 
repeated for each data point, and the final data value will be -100000.00, 
figure A28. The following cards after the final data value will be the 
same as for ADAS - Radar Combination Data. 
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ADAS - Radar Combination Data 
System Card Data 

Radar Card Data 

Prints ADAS and Radar Data 
before merging 

Prints data written onto 

tape 16 (NFF) 

Produces Calcomp plot of data 
on tape 16 (Calcomp plot #1) 
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Termination of card by a $ sign. 


* This card is of the NAMELIST format. 
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** PHDG = 1 for a heading. 
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(R50, CEP, confidence limits) 


Y-axis length (inches) for plot no. 2 


Starting value of Y-axis for plot no. 


Starting value for X-axis is the 


same for all 3 plots. 
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latitude and longitude data 


are in degrees. 
System Card Data 
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lope of CEP Specification line for 


s 
Calcomp plot 





Termination of card by a $ sign. 


SLOPE=XX.XXXX (nm/unit time of plot absissa 


This card is of the NAMELIST format. 


NOTE: 





* Input Options are the same as for Program 
NAVAN Data Card #1. 
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Elapsed time (seconds) 
Latitude error (degrees) 
Longitude error (degrees) 


Flight number 
Figure A23 
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Final data value 


-100000.00 
Figure A24 
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INTRODUCTION 


The Circular Error of Probability Program (CEPLOT) calculates the 
CEP and plots the data. The input is navigation system end point position 
error data normalized to one hour, based on a linear error rate. Data are 
input in card form. 


OVERVIEW 


CEPLOT consists of the main program and subprograms. It reads the 
data, computes the error, normalizes the data, computes the CEP, and 
produces a tabular listing and Calcomp plot tape for plotting. 


A flowchart is shown in figure Bl- 


The method used to calculate the CEP is the same as is used in NAVAN 
(Appendix A), which is based on the Air Standard 53/11.B, 15 August 1968, 
The Specification and Evaluation of the Accuracy of Inertial Navigation 
System. 


PREPARATION FOR USE 


The program deck of CEPLOT will be permanently stored at the AFFTC 
Systems Engineering Branch. 


USER'S GUIDE 


The deck developed for running CEPLOT is composed of three parts. 
The first part is the job control cards which are used to attach the 
Calcomp plot tapes and execute the program. Examples of these control 
cards are shown in figure B2. 


The second portion of the program is the main program and 
subroutines. 


The third portion of the program consists of two parts: (1) Tne 
parameter cards which control the program and plotting, and (2) The data 
cards. These data cards are shown in figure B3. 


The output data from CEPLOT consists of a computer output listing, 
figure B4, and a Calcomp plot, figure B5. 


INPUT CARDS FOR CEPLOT 


CEPLOT data card #1, figure B6, contains the number of data points to 
be calculated and plotted, and the heading to be placed on the Calcomp 
plot. 


CEPLOT data card #2, figure B7, contains the type of INS alignment, 
either gyrocompass (GC), best azimuth heading (BATH), or standard (STD), 
and the format of the input data. The data format may be "RAW" which 
means that the delta latitude, delta longitude, flight number, and time r 
will be input. Radial error will be calculated and the data normalized. 
In the other format, the data are input as delta latitude, delta longitude, 
and radial error (PMT Left Blank). j 


CEPLOT data card #3, figure B8, the first data card, if not "RAW", 
will contain delta latitude, delta longitude, and radial error which has 


a 





been normalized. One card is used for each data point. If "RAW" data, 
the card, figure B9, will have flight number, delta latitude, delta 
longitude, and time in hours and minutes. This data will not have been 
normalized. One card is required for each data point. 
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Check Cases 








The following check case for program NAVAN uses data from an 
onboard navigation system and space positioning (RADAR) tracking data. 
The ADAS and RADAR data are first merged and placed on the New History 
Plight File, then another set of data is added to this file and the 
merged data is analyzed and the results are plotted. All intermediate 
data was printed for help in diagnosis of program errors. 
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The check case for program CEPLOT is a set of random end 
errors (generated using a random number routine) and illustrat 


program operation including plotting of the data. 
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azimuth 
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GLOSSARY 


elastic deformation of a gyro's gimbals that result 
in gyro drift proportional to the g? value of 
vibration. 


horizontal direction or bearing. 


azimuth measured from 0° at the north or south 
reference direction clockwise or counterclockwise 
through 90° or 180. 


a statistical significance test based on frequency 
of occurance. 


a measure of the accuracy with which an aircraft 
can be guided; the radius of the circle ata 
specific distance in which 50 percent of navigation 
errors fall; also called circular error probable, 
and circle of probable error. 


The position of the accelerometer pendulum is off 
null by a small angle in order to create the error 
Signal necessary to rebalance the accelerometer. 
This angle causes cross-coupling of the component 
of acceleration which is normal to the 
accelerometer. 


a measure of central position. The geometric mean 
of n quantities equals the nth root of the product 
of the quantities. 


horizontal direction in which an aircraft is 
pointed, expressed as angular distance from a 
reference direction. 


a system in which the (vector) momentum of a 
particle is conserved in the absence of external 
forces. Thus, only in an inertial system can 
Newton's Laws of Motion be appropriately applied. 


Terrestrial latitude is angular distance from the 
Eqyator, measured northward or southward through 
90°. 


Terrestrial longitude is the arc of a parallel, or 
the angle at the pole, between the prime meridian 
and the meridian of a point on the Earth, measured 
eastward or westward from the prime meridian 
through 180°. 


magnetic pole 


Proportional bias 


root mean square 
error (RMS) 


Student's t 
distribution test 


vibropendulous 








either of the two places on the sux face of the 
Earth where the magnetic dip is 90°, that in the 
Nogthern Hemisphere (at, approximately, latitude, 
73°8°N, longitude, 101°W in 1955) being designated 
North Magnetic Pole, and that in the Southegn 
Hemisphere (at, approximately, latitude, 68 S, 


longitude, 144°E in 1955) being designated South 
Magnetic Pole. 


the bias required to prevent gyro drift that 
results from acceleration perpendicular to the 
output axis at the gyro. 


in statistics, the square root of the arithmetic 
mean of the squares of the deviations of the 
various items from the arithmetic mean of the 
whole, also termed standard deviation. 


a statistical method to test the hypothesis that 
the mean of the sample is consistant with being 
equal to the mean, yu, of the assumed population. 


the accelerometer error produced when the 
pendulous mass in an accelerometer is vibrated. 
This error is proportional to the vibropendulous 
coefficient times this vibration in feet per 
second per second. 
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LIST OF ABBREVIATIONS AND SYMBOLS 


Definition 

acceleration 

centripetal acceleration 
distance 

East 

acceleration of gravity at mean sea level 
geometric nean . 
(subscript) test, i=1,m 
confidence limits 

number of tests 

North 

a variable 

sample mean of q 

radius of the Earth 

radial error = Vv x? + y? 
GM/RMS 

pth percentile of radial error 
South 

time 

t in Student's t test 
Technical Order 

velocity 


Vest 


‘latitude error 


longitude error 
vertical 


pth percentile point of a zero mean 
normal distribution 






Unit 


ft per sec? 


ft per sec? 


NM or ft 


32.3 ft per sec? 


dimensionless 
dimensionless 


dimensionless 


dimensionless 


dimensionless 


dimensionless 


2.09 x 10” ft 
NM 
dimensionless 


dimensionless 


ft per sec 
NM 
NM 


dimensionless 








Definition 


(l-a)% confidence intervals 
azimuth drift rate 
acceleration error 

distance error 

level gyro drift rate 
velocity error inside the Schuler loop 
velocity error outside the Schuler loop 
azimuth angle error 

pitch angle error 

roll angle error 

pitch angle 

a/r 

mean or expected value of q 
sample sigma 

sample variance of q 
azimuth angle 

variance of q 

sum of qs for i=i1,m 
latitude 

roll angle 

chi 

gyro torquing rate 


Earth rotation rate 


radians per unit time 





Unit 
dimensionless 
deg per hr 

ft per sec? 
NM or ft 

deg per hr 


ft per sec 


ft per sec 
deg 

deg 

deg 

deg 

1/sec? 
dimensionless 
dimensionless 
dimensionless , 
deq 
dimensionless 
dimensionless 
deg 

deg 
dimensionless 
deg per hr 


15.04 deg per hr 


rad per hr 





